ASTROPHYSICAL JOURNAL, IN PRESS 

Preprint typeset using lAT^^X style cmulateapj v. 04/03/99 



INTERGALACTIC DUST AND OBSERVATIONS OF TYPE lA SUPERNOVAE 

Anthony Aguirre 

Department of Astronomy, Harvard University 
60 Garden Street, Cambridge, MA 02138, USA 
email: aaguirre@cfa.harvard.edu 
Astrophysical Journal, in press 



a^ 
a^ 

1—5 

(N 

m 
> 

m 
o 

6 • 



ABSTRACT 

Estimates of the cosmic star formation rate and of cluster metallicities independently imply that at 
z <^ 0.5 the gas in the universe has substantial average metallicity: 1/10 <, Z/Zq <5 1/3 for flgas — 0.05. 
This metal density probably cannot be contained in known solar-metallicity galaxies of density parameter 
ri* « 0.004, implying significant enrichment of the intergalactic medium (IGM) by ejection of metals 
and dust from galaxies via winds, in mergers or in dust efflux driven by radiation pressure. Galaxies 
have a dust/metal ratio of 0.5 in their interstellar media, but some fraction (1 — /) > of this must 
be destroyed in the IGM or during the ejection process. Assuming the Draine & Lee dust model and 
preferential destruction of small grains (as destruction by sputtering would provide), I calculate the 
reddening and extinction of a uniform cosmological dust component in terms of / and the minimum 
grain size amm ■ Very small grains provide most of the reddening but less than half of the opacity for 
optical extinction. For / ^ 0.3 and amin ^ 0.1 /im, the intergalactic dust would be too grey to have 
been detected by its reddening, yet dense enough to be cosmologically important: it could account for 
the recently observed type la supernova dimming at z ^ 0.5 without cosmic acceleration. It would also 
have implications for galaxy counts and evolutionary studies, and would contribute significantly to the 
cosmic infrared background (GIB). The importance of grey intergalactic dust of the described type can 
be tested by observations of 2 = 0.5 supernovae in (rest) i?— band or longer wavelengths and by the 
fluxes of a large sample of supernovae at z > 1. 

Subject headings: cosmology: observations, distance scale - intergalactic medium - dust, extinction 



1. INTRODUCTION 

The importance of dust extinction in the Galaxy has 
been recognized since early in this century when star- 
counting surveys revealed absorption of optical light by 
'dark clouds' (Barnard 1919). It is fortunate that extinc- 
tion correlates relatively well with reddening in the Galaxy, 
because it is difficult enough to accurately measure either 
the distance to a typical astronomical source or its intrin- 
sic luminosity - let alone both. But knowing the intrinsic 
color {B — V)i (using an unobscured line of sight) along 
with the observed B — V color and a reddening-extinction 
law Av = Rv[{B -V)-{B- V)i] = RvE{B - V), one 
can correctly determine the distance of an object from its 
distance modulus, or vice versa. The value of Ry varies 
markedly within the Milky Way {i <, Rv <, 6; Mathis 
1990) and among different galaxies (1.5 ^ Rv <, 7.2; Falco 
et al. 1999), but a value of Ry — 3.2 is useful for many 
estimates of Galactic extinction. 

Radiation from extragalactic objects is subject to ex- 
tinction by dust both inside and outside of the Galaxy. 
However, while extragalactic dust has received attention, 
our knowledge of its amount and properties is rather lim- 
ited, because methods useful for estimating dust density 
in galaxies have proved less effective when extended to in- 
tergalactic space. For example, a number of groups have 
attempted to measure extinction by dust in clusters us- 
ing background-object counting, and several claims of in- 
tracluster dust (Bogart & Wagoner 1973; Boyle, Fong & 
Shanks 1988; Romani & Maoz 1992) and extragalactic 
dust cloud detections (Wsozlek et al. 1988 and references 



therein) have been made, but even now these claims re- 
main controversial (see Maoz 1995). The strongest pro- 
posed limits on a diffuse distribution of intergalactic dust 
with a Galactic reddening law have come from studies of 
the redshift evolution of the mean quasar spectral index 
(e.g. Wright 1981; Wright & Malkan 1988; Cheng, Gaskell 
& Koratkar 1991). These studies limit uniform dust of 
constant comoving density to have Av{z = 1) < 0.05 mag 
(from Wright & Malkan 1988), and are most sensitive to 
dust at z > 1. 

While our knowledge of it is poor, intergalactic dust 
could have great cosmological importance, as it could af- 
fect results concerning the cosmic microwave (CMB) and 
cosmic infrared (CIB) backgrounds, galaxy and quasar 
numbers at high z, galaxy evolution, large-scale structure, 
etc. This paper discusses intergalactic dust chiefly in the 
context of its importance in measurements of the cosmo- 
logical deceleration parameter - a subject discussed nu- 
merous times, first by Eigenson (1949) and most recently 
in Aguirre 1999 (A99). 

Conditions in the diffuse intergalactic medium (IGM) 
strongly disfavor dust formation, so whatever intergalac- 
tic dust exists is probably either the remnant of an early 
Population III epoch, or is formed in galaxies and removed 
by some mechanism. (The remaining possibility, that a 
substantial dust-forming population of extragalactic stars 
exists, is not considered here.) Previous investigations of 
intergalactic dust have almost invariably assumed that it 
has properties similar to that of Galactic dust; but this 
assumption is not well justified. Even among galaxies, Ry 
varies by a factor of four, and intergalactic dust may have 



2 



INTERGALACTIC DUST AND TYPE lA SUPERNOVAE 



creation, destruction, and selection mechanisms quite dif- 
ferent from dust in galaxies. As argued in A99, radiation 
pressure ejects grains with high opacity and a broad opac- 
ity curve more efficiently than other grain types. In §^ I 
discuss results suggesting that small grains are preferen- 
tially destroyed by sputtering, both in the halos of galaxies 
(during the ejection process) and perhaps in the IGM. 

The large grains, while having significant mass, give very 
small E{B — V) reddening, and actually have higher visual 
opacity (per unit mass) than dust which includes small 
grains. Intergalactic dust of this type would not have been 
detected by quasar reddening surveys. Current data, de- 
scribed in §^ suggests that the universe has been enriched 
to ^ 1/10 solar metallicity before z ^ 0.5. Section [s] shows 
that if a significant fraction of this metal is locked in dust 
that is distributed fairly uniformly, the dust extinction to 
z ~ 0.5 can explain the observed progressive dimming of 
type la supernovae (Riess et al. 1998; Perlmutter et al. 
1999 [P99]) without cosmic acceleration. 

The last section outlines ways in which the influence of 
the type of intergalactic dust described here can be tested, 
probably most cleanly by future supernova observations. 

2. COSMIC METALLICITY 

Although still very incomplete, our understanding of the 
evolution of cosmic metallicity has improved dramatically 
during the last few years. Recent surveys (see e.g. Madau 
1999 for a review) of galaxies in the Hubble Deep Field 
show that the comoving star formation rate (SFR) rises 
with redshift by a factor of ten to z ~ 1 — 1.5, past which 
it either declines or levels off. Using the observed star 
formation rate and an assumed ratio of metal formation 
to star formation of 1/42 (Madau et al. 1996), the total 
cosmic metallicity (neglecting any population III contri- 
bution) may be estimated by integrating over time the 
metal formation rate. Figure shows the result of this 
integration starting at z = 10, using curves from Madau 
(1999) and Steidel et al. (1999). These demonstrate that 
whether or not the SFR declines for z ^ 1.5, the cur- 
rent metal density is ~ (1-5 — 2)hQ^ x ID""*. The es- 
timates shown for flz{z) agree with other results in the 
literature using similar methods: Pettini (1999) estimates 
nz(z~ 2.5) - 6 X 10"^ Gen & Ostriker (1999b) calculate 
riziz ~ 2.5) - 2 X 10-^ and flziz ~ 0.5) ~ 1.1 x IQ-*. 

Very interestingly, these estimates coincide with the 'fos- 
sil evidence' presented by Renzini (1997; 1998). He argues 
that clusters are essentially closed systems which contain 
all of the metal produced by their stellar populations. 
Stars of approximately solar metallicity comprise only a 

3/2 

fraction Td ~ 0.09/ig5 of the total cluster gas mass, yet 
the remaining intracluster gas has w 1/3 solar metallic- 
ity. The associated metal production per unit of stellar 
mass can be written Mz ~ M^[l + 3.15hQ^^'^]ZQ (Renzini 
1997). Unless stars in clusters produce metals much more 
efRciently than those in field galaxies, this figure should ap- 
ply to the « 0.004 (Fukugita, Hogan & Peebles 1996) 
of stars in the universe, giving flz « 3.1 x 10"'*. More- 
over, if the cosmic gas density is flgas ~ 0.05, then the 
star formation efficiency of clusters well represents that of 
the universe, Tigm = rP* ~ Tc;, further indicating that 



clusters are a fair sample. 

Two obvious ratios link these numbers to estimates of 
intergalactic dust density: the fraction Fj of all metal re- 
siding in the IGM, and the fraction dm of the intergalactic 
metal locked in dust. At high z, the primary data bear- 
ing upon these ratios are observations of metal lines in 
quasar absorption spectra. Pettini et al. (1997) have mea- 
sured the depletions of Gr and Zn in damped Lyman-a 
systems which Pei, Fall & Hauser (1998) interpret as giv- 
ing a dust/metal ratio of « 0.5 for all z ^ 3, similar to 
that of known galaxies. This suggests that even at high 
z, source regions for the cosmic metallicity have dm ~ 0.5. 
Whether this applies to metal as it reaches the IGM de- 
pends on the degree to which the ejection process destroys 
dust, or selectively ejects either gaseous metal or dust (to 
be discussed in §2). 

Observations of lower column density Ly-a absorbers 
give information about the gaseous metals in relatively 
cool regions of the IGM. Under some (important) assump- 
tions about ionization corrections, these line strengths 
indicates that areas with column density N{HI) « 
-|^q14-15 jjave ~ 0.3% solar metallicity, giving 

rj'|(z - 3) - 3 X 10-6 in gaseous metal if this enrich- 
ment is uniform. Direct integration (effectively treating 
N{HI) <; lO^^'cm^^ regions as pristine) gives fi'^(z ^ 
3) > 3 X 10-"^ (Songaila 1997). In hght of the higher num- 
bers derived from the SFR, the assumption that most of 
the cosmic metallicity resides in the Ly-a forest therefore 
leads to a 'missing metals' problem, as noted by Pettini 
(1999) and Renzini (1999), who argue that this discrep- 
ancy may indicate that most metals are located in hot 
halos around galaxies, in proto-clusters, or in a phase of 
the IGM hotter than that which the Ly-a forest studies 
probe. Gen & Ostriker (1999a) reach the same conclusion 
on separate grounds, arguing on the basis of both simple 
physical arguments and their numerical simulations that 
the bulk of universal baryons at low z should be hot. 

The idea that the IGM is fairly metal-rich gains more 
support from arguments that hold at low redshift. If en- 
riched to solar metallicity, the current mass of stars and 
gas in known galaxies is sufficient to contain the metals ex- 
pected at z ^ 2.5 (if most galaxies had already formed by 
then), but this does not hold for later epochs: as also ar- 
gued by Renzini (1998), the cosmic metallicity at z = de- 
rived from cluster enrichment or by integration of the SFR 
cannot be stored in currently observed galaxies, which can 
hold at most 
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The remaining « 50 — 75% must exist in the intergalactic 
gas, in extended halos, or hidden in undetected galaxies 
(which seems unlikely). This is exactly what one would 
expect on the basis of the metal distribution in clusters: 
Renzini (1997) argues that while intracluster gas contains 
about three times as much metal as the cluster galaxies, 
there seems to be no compelling reason to expect that 
metal ejection is much more efficient in clusters than in 
the field^ , so most metals created from field galaxies should 
also lie outside them. 

^Even if ram-pressure stripping and/or mergers made clusters more efficient, this would predict that field galaxies would be about four 
times as metal-rich as cluster galaxies, contrary to observations. Also, clusters with higher velocity dispersion should have higher metallicity; 
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To avoid the conclusion that the universe has substan- 
tial intergalactic nietallicity, it therefore seems that one 
would have to argue both that cluster galaxies eject met- 
als more efficiently than field galaxies and have a different 
IMF, and also accept that estimates of the SFR and/or 
the density of metal in galaxies are incorrect by at least a 
factor of two. 

In summary, the density of intergalactic dust can be es- 
timated at 2 ~ 0.5 as 

~ 0.5) =Fjxd^xxx 10-^ (2) 

with likely values of 1.5 ^ x <; 3 and — 0.5. The argu- 
ment that stars and gas in known galaxies cannot contain 
these metals gives 0.5 ^ Fj ^ 0.75; the larger number also 
corresponds to the value derived by assuming that field 

galaxies eject metal as cluster galaxies do. The resulting 
estimate does not take into account grain destruction. 

3. REMOVAL OF DUST AND METAL FROM GALAXIES 

The conclusion that the IGM is fairly metal rich implies 
that metal can be efficiently removed from the galaxies in 
which it forms. Several ways in which galaxies can eject 
dust and metallic gas have been studied in the context of 
clusters (for exactly the same reason), but it rather unclear 
which mechanism is dominant. The density of intracluster 
gas is high enough that ram-pressure stripping of galactic 
gas may be efficient (e.g. Fukumoto & Ikeuchi 1996; Gunn 
& Gott 1972) in clusters, but it would be much less effec- 
tive for a galaxy in the general IGM. The removal of gas by 
supernova-driven winds has been widely discussed, and de- 
tailed simulations have been performed investigating this 
effect in large starburst galaxies (Suchkov et al. 1994) and 
in dwarves (Mac Low & Ferrara 1998). It is widely thought 
that galactic winds also remove gas from ellipticals (e.g. 
David, Forman & Jones 1990). Gnedin (1998) has claimed 
that mergers provide the dominant metal removal mech- 
anism, at least for z ^ 4. Finally, dust removal (without 
metallic gas) by radiation pressure can be fairly efficient 
even for present-day spirals; starburst galaxies with higher 
luminosities would be correspondingly more effective. In 
this section I will concentrate on radiation pressure as the 
dominant dust expulsion mechanism because the effects 
on the dust have been investigated most carefully in that 
scenario; but I will discuss the other mechanisms briefly. 

The dynamical evolution of a dust grain in a spiral 
galaxy is governed primarily by the radiation pressure, 
gravity, the viscous gas drag and the magnetic Lorentz 
force. Starting with Chiao & Wickramasinghe (1972), sev- 
eral groups have studied these forces acting on grains in 
model galaxies with some assumptions about the mass, 
gas, and luminosity distributions. Confirming the results 
of Chiao & Wickramasinghe, Ferrara et al. (1990) find 
that graphite grains of most sizes can escape most spiral 
galaxies, and that silicate grains are marginally confined 
(though silicate grains with a ^ 0.05 /xm can escape high 
luminosity spirals.) 

Calculating the grain dynamics for two specific galactic 
models (of the Milky way and NGC 3198) in more de- 
tail, Ferrara et al. (1991) find that in the Milky Way, 
silicate grains of radius a = 0.1 /xm have typical speeds 



of ^ 200 — 600 km s^^ and reach halo radii ^ 100 kpc in 
^ 150 — 500 Myr. Graphite grains of a = 0.05 /im can 
move approximately twice as fast. Another investigation, 
by Shustov & Vibe (1995), gives similar results. They find 
that grains of size 0.07 /im < a < 0.2 /im are ejected most 
effectively. Silicate (graphite) grains of 0.1 /im starting 
1 kpc above the galactic disk attain speeds of 1000 km s^^ 
(2000kms-i) and reach 100 kpc in 100 (40) Myr. Shustov 
& Vibe stress that only dust starting somewhat above the 
disk can escape, but this does not imply that dust ejec- 
tion is inefficient; in their model, the dust expulsion rate 
is up to O.4M0yr~^, which in a Hubble time exceeds the 
entire metal content of the Galaxy and gives ^i^^™ of or- 
der lO"'^, even assuming a constant SFR. Most recently, 
Davies et al. (1998) have performed numerical calculations 
of dust removal by radiation pressure, taking into account 
the opacity of the disk. Disk opacity reduces the radi- 
ation pressure at high galactic latitude, so Davies et al. 
find dust expulsion less efficient than the previous stud- 
ies.^ Nevertheless, Davies et al. predict the removal of 
(at least) 0.1 /zm graphite grains from their model galaxy 
if there is fairly low disk opacity. Smaller grains are ex- 
pelled much less efficiently. Notably, Davies et al. still 
estimate that up to 90% of the dust formed in spirals may 
be ejected, and even calculate an estimate of 1 mag of 
intergalactic extinction across a Hubble distance. 

All studies of radiation pressure driven dust removal 
have noted the importance of magnetic fields on grain dy- 
namics but have (effectively) neglected them in their cal- 
culations, with the following justifications given: 

• Magnetic field lines are only sometimes parallel to 
the disk, and only on large scales (all studies). 

• Grains are charged only part of the time (Ferrara et 
al. 1991; Davies et al. 1998). 

• Radiation pressure can enhance Parker (1972) in- 
stabilities that can lead to open field configurations 
with lines perpendicular to the disk. (Chiao & Wick- 
ramasinghe 1972; Ferrara et al. 1991). 

An additional justification is empirical: dust is actually 
observed outside the disks of galaxies (Howk & Savage 
1999; Alton, Davies & Bianchi 1999; Davies et al. 1998 
and references therein; Ferrara et al. 1999 and references 
therein). However the dust escapes the disk, its presence 
proves that while magnetic fields are potentially important 
and currently impossible to model in detail, they cannot 
be perfectly effective at dust confinement. On the other 
hand, this does not prove that dust can fully decouple from 
the gas. 

While dust expelled by radiation pressure could not 
carry a significant gas mass with it, other metal expulsion 
mechanisms probably remove dust along with gas. Alton 
et. al. (1999) have presented observations of dust outfiows 
in three nearby starburst galaxies, concluding that the 'su- 
perwinds' driving these outflows can impart near-escape 
velocity on the dust, and that up to 10% of the dust mass 
of these galaxies could be lost in the observed outflows 



Renzini finds no evidence for this. 

^Of course a higher assumed intrinsic luminosity of the galaxy would cancel this effect; Davies et al. only investigate one mass-to-light ratio. 
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alone. These results lend observational support to the nu- 
merical simulations of Suchkov et al. (1994) which pre- 
dicted such outflows, and further demonstrate that dust 
can escape along with metallic gas. Lehnert & Heckman 
(1996) have estimated the efficiency of metal removal by 
winds in starburst galaxies using a large sample, and find 
that galaxies could enrich the IGM to O^™ ^ 5 x 10^^. 
This figure assumes a constant SFR and the authors esti- 
mate that the enrichment is likely to be ten times higher 
with a more realistic SFR. 

Supernova-driven winds might also eject dust from ellip- 
ticals. Vereshchagin, Smirnov & Tutukov (1989) estimate 
that the ratio of galactic wind force to gravitational force 
on a grain of radius a is 
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Vw/a, 



where Vw is the wind velocity and a is the specific mass 
ejection rate for the wind in the galaxy. For a ~ 10~^^ s^^ 
(applicable for the present epoch; Vereshchagin et al. 
1989; David et al. 1990) and relatively slow winds {Vw ^ 
10 — 60kms~^; Vereshchagin et al. 1989), only very small 
grains escape; but for starburst ellipticals, David et al. find 



sharp difference arises because for sputtering at fixed gas 
temperature and grain velocity the destruction timescale 
{l/a){da/dt) oc 1/a (Draine & Salpeter 1979a), and be- 
cause small grains are more affected by gas drag yet less 
propelled by radiation pressure, hence move more slowly 
through the halo.'^ Shustov & Vibe conclude from their 
calculations that the grains escaping intact from galaxies 
will have sizes 0.03 fim <, a <^ 0.2 fun for graphite par- 
ticles and 0.07 /im <^ a < 0.2 /im for silicate particles.'* 
These conclusions depend on assumptions about the den- 
sity, temperature, and extent of the galactic halos, but the 
fact that both groups obtain similar results suggests that 
the minimal grain size surviving expulsion is probably of 
order amin 0.05 ^m. 

The efficiency of dust destruction in other metal re- 
moval processes has not been calculated in detail and is 
difficult to estimate. Dust driven out by winds would be 
vulnerable to sputtering by the halo gas as well as by the 
faster moving wind, though it will be somewhat shielded 
if embedded in cool clumps of gas. To estimate the effect 
of the wind, let us assume a mass loss rate M (in solar 
masses/yr) due to a wind leaving the galaxy radially with 
velocity Vw- The effect of this wind would be similar to 



the effect of a hot gas of temperature T] 



and (proton) number density Hp 3M /IGirR'^Vw'rnp. For 
125 km s"^ < V ^ 4000 km s^^ and B ^ 10 kpc, this gives 

10^ K < T < 10^ K and 



1.9 X IQ-^M 



Using Draine & Salpeter's (1979a) 
graphite in this temperature range, 
a lifetime of 
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rates of at least 5x (lO-^^-lO^i^) s"* for the first lO^yr of 
starburst activity, implying that the winds dominate grav- 
ity in grain dynamics for grains up to ~ 0.05 — 2.5 ^m, even 
for winds too slow to escape the galaxy themselves. Es- 
sentially the same argument would apply to grains subject 
to winds in spirals, so whether dust lies in relatively cool 
outflows or is exposed to the wind itself, it is difficult to 
see how it could avoid being driven out with the metallic 
gas. 

Gnedin (1998) has performed high resolution cosmolog- 
ical simulations that suggest that mergers are the domi- 
nant metal removal mechanism, at least at high (z ^ 4) 
redshift. This mechanism would eject gas with roughly the 
same dust/gas ratio as the source galaxy. 

This brief survey of metal ejection mechanisms suggests 
that it is difficult to efficiently remove metallic gas from 
galaxies without also removing dust (although the con- 
verse of this would not be true if radiation pressure is the 
dominant ejection mechanism). Rough estimates of the 
ejection efficiencies show that metal ejection rates suffi- 
cient to account for the enrichment of the IGM estimated 
in §0 are reasonable. While dust probably accompanies 
gas as it leaves galaxies (or leaves by itself), studies of 
clusters show that intracluster gas is not dust-rich. The 
next section addresses the probable cause of this disagree- 
ment: grain destruction during the ejection process and in 
the IGM. 

4. DESTRUCTION OF SMALL GRAINS 

A key result of the investigations by both Ferrara et 
al. (1991) and Shustov & Vibe (1995) is that grain de- 
struction due to sputtering by hot halo gas is relatively 
insignificant for grains of a ~ 0.1 /im but very effective 
for grains with a ^ 0.01 ^m. While 0.1 /xm grains lose 
only ~ 0.005 /zm in radius, the small grains are com- 
pletely destroyed on a timescale of ^ 500 Myr. The 

^Size effects can be even stronger; Draine & Salpeter (1979b) find that the most efficient dust destruction 
medium', is ~ 500 X more effective in 0.01 fim grains than for a = 0.1 ^m. 

*These numbers are rather approximate because the authors computed results only for six grain radii. 
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sputtering rate for 
this corresponds to 



Vw 



125 km s- 



0.01 /im 



This is comparable to the ejection timescale, so this sput- 
tering could be important but is unlikely to completely 
destroy the dust. 

Grain-grain collisions provide another important dust 
destruction mechanism in galaxies and might be impor- 
tant in the early stages of the ejection process. For ex- 
ample, grain-grain collisions in supernova shocks can effi- 
ciently shatter large grains into smaller ones (e.g. Jones, 
Tielens & Hollenbach 1996), so if supernova blowout re- 
moves dust, there is a danger that shocks from the same 
supernovae might shatter the large grains before the dust 
is expelled. Shocks may also play an important role in 
mergers. On the other hand, it is unclear whether the 
dust observed in the ISM is representative of dust which 
has just formed, or already been shock-processed, or some 
steady-state between the two. Specifically, there is evi- 
dence for the formation of large grains in novae (Shore et 
al. 1994), and possibly in supernovae (see Wooden 1997 
and Pun et al. 1995), and grains are presumably larger in 
molecular clouds where high values of Ry are measured. It 
may be, then, that pre-shock grains tend to be somewhat 

sputtering in the 'inter-cloud 



AGUIRRE 



5 



larger, and the MRN distribution is more characteristic of 
grains after significant shattering has occurred. The as- 
sumption of this paper is that dust leaving its progenitor 
galaxy will have an grain size distribution characteristic of 
dust in the ISM. In the absence of significant shattering, 
this is probably conservative, since a significant fraction of 
dust is contained in dense clouds with high Ry (e.g. Kim, 
Martin & Hendry 1994). 

4.1. Dust Destruction in the IGM 

Rather little is known about the destruction of dust 
in the IGM. Schmidt (1974) estimates that soft cosmic 
rays would provide the most efficient destruction, but can- 
not determine whether or not the destruction time would 
exceed the Hubble time; moreover, Draine & Salpeter 
(1979b) find that cosmic rays are unimportant dust de- 
stroyers in the Galaxy (where they should be at least as ef- 
fective as in the IGM). The hot gas component of the IGM, 
however, could sputter grains effectively, even at low den- 
sity. Using again Draine & Salpeter's (1979a) estimate^ , 
the lifetime can be written 



and details of the dominant mechanism of metal ejection 
for galaxies are uncertain. Thus the arguments of this 
section are intended merely to make plausible the chief 
assumption of this paper, which is that grains of size 
a <, 0.05 — 0.1 /im are removed (either by destruction or by 
failure to escape their progenitor galaxies) from the grain- 
size distribution characterizing dust outside of galaxies, 
whereas larger grains are not. 

5. DENSITY OF SURVIVING INTERGALACTIC DUST 

The estimate of the density of intergalactic dust in sec- 
tion §0 did not take into account dust destruction or the 
preferential expulsion of dust. Lets us assume that a mass 
fraction (1 — /esc) of dust is destroyed as it leaves the 
disk and/or traverses the halo, and that a further fraction 
(1 — figm) is destroyed in the IGM after the dust escapes 
the halos but before z ~ 0.5. There are three general 
scenarios indicated by the dust ejection and destruction 
mechanisms outline above: 

1. Dust and gas leave together, with the galactic 



(3.5-9) X 10^ 



0.01 /im 



.... 

where ^Igas signifies the hot gas density in critical units 
and ^ is a clumping factor. The Hubble time (for = 1) 
is H~^{z) = 1.6 X 10^°/ig5^(l + z)^^/^yr, suggesting the 
efficient destruction of grains for which 



dust/metals ratio of « 0.5 
this survives, so that dm ^ 



A fraction fescfrgm of 
This sce- 



g EE 0.14 



0.01 /im 



^gas 

005 



/l65^(l + ^)'/' » 1. (5) 



^ -^f esc figm- 

nario predicts a high enrichment of the IGM near 
galaxies and perhaps a substantial density of metal 
in galactic halos. 

2. Dust, driven by radiation pressure, decouples from 
the gas either in the disk or in the inner halo, but 
is partially destroyed. The gaseous metal (both de- 
stroyed dust and metal which escapes the disk but 
then decouples from dust) could (a) remain in the 
halo or could (b) return to the disk to form more 
dust, repeating the process. In the former case 
galactic halos may be highly enriched and dm ~ 
Q.^f esc figm] in the latter (unlikely) case galaxies 
should be very deficient in metals which are easily 
incorporated into dust, and dm figm is possible. 

3. Gaseous metals leave disks but dust remains {dm ^ 
1). While unlikely, this would lead to a highly en- 
riched IGM and/or halo gas component but little in- 
tergalactic obscuration (like the case f esc figm ^1). 
Disks would be heavily enriched with elements that 
do form dust. 

Assuming that 0.5 < Fj <^ 0.75, equation |2| gives 

j(z ~ 0.5) ~ (4 - 11) X 10-^ fescf^gm 



Qtgm 
"dust*- 
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The clumping factor (i.e. the overdensity felt by a 'typi- 
cal' grain) is quite uncertain, but the simulations of Gen 
& Ostriker (1999b), which numerically track the distri- 
bution of metallicity, indicate that at z 0.5, the mean 
universal metallicity approaches the metallicity of ^ ^ 100 
regions. Regions of much higher overdensity do not have 
much higher metallicity and hence cannot contain most 
of the metals - for example, S ^ 1000 regions have only 
about twice the metallicity, so dense 'subregions' can con- 
tain only about 20% of the metals in 5 ~ 100 regions. If 
a 'typical' grain experiences 5 ~ 100, then Q ^ 26 for 
0.01 /im grains and Q ~ 2.6 for 0.1 fiui grains. This is 
suggestive (but only suggestive) that sputtering by hot in- 
tergalactic gas might provide yet another mechanism by 
which grains of a ^ 0.1 fim might be selectively destroyed. 

Finally, note that the low mean dust density in the IGM 
and in extended galaxy halos would strongly suppress the 
grain-grain collisions thought to shatter large grains into 
small ones in the galaxy^; since dust formation is also in- 
efficient in the IGM there is probably no source of new 
small grains outside of galaxies. 

The efficiency of dust destruction depends in a rather 
complicated way on the environment; moreover the type 

^More recent calculations by Tielens et al. give similar sputtering rates for carbon at T > 10^ K, while their rates are somewhat higher for 

silicates and somewhat lower in both materials for 10^ K < T < lO'^ K. 

^Equation ^ gives the dust optical depth through the halo of a galaxy. The 'optical depth' to an emerging grain would be of similar 
magnitude, so grain-grain collisions are probably unimportant unless high-z galaxies are all heavily obscured by dust in their halos. 

^Sputtering will also change the upper grain-size cutoff, in effect shifting the whole distribution toward smaller radii. The neglect of this 

fajt may be justified by the excess of large grains over the MRN prediction indicated by estimates of the actual grain-size distribution (see 
ih, and because sputtering may be more selective in destroying small grains than the 1/a behavior would imply. 



for scenarios 1 and 2a. The dust density would be higher 
by a factor of up to 2//esc for scenario 2b, and very 
small for 3. 

Since small grains are preferentially destroyed, but prob- 
ably cannot be created in halos and in the IGM, fescfigm 
effectively determines the minimal grain size amin-"^ The 
next section discusses the dependence of dust properties 
on ttmin, and gives corresponding values of fescfigm- 
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6. PROPERTIES OF THE DUST 

The absence of a <^ 0.1 /im grains would have impor- 
tant impUcations for the properties of intergalactic dust. 
The most commonly used model for Galactic dust is the 
two component Draine & Lee (1984; DL) model consisting 
of silicate and graphite spheres with a distribution in ra- 
dius (as proposed by Mathis, Rumpl and Nordsieck 1977; 
MRN) oiN{a)da oc a"3■^ 0.005 ^im <a< 0.25 /xm. After 
synthesizing dielectric functions for both graphite and 'as- 
tronomical silicate' and assuming a silicate/graphite mass 
ratio ~ 1, DL demonstrated that the resulting model fits 
both the observed opacity and polarization over a wide 
wavelength range (0.1 /im <^ A < 1000 /xm), most no- 
tably fitting the observed features at 0.2175 /xm and 10 ^m. 
This paper employs the DL model not because it is most 
likely to be correct, but because there is little agreement 
as to what the correct grain model might be. The DL 
model is widely used and familiar, and hopefully (but by 
no means certainly) captures the essential features of the 
dust. Other models are discussed briefly below. 

An interesting aspect of the MRN distribution is that 
while the geometrical cross section (ex a^) is dominated 
by small-radius grains, the mass (cx a^) is dominated by 
grains of large radii. Thus, removing the very small grains 
can affect the opacity curve dramatically, without radically 
changing the total dust mass. 

Figure ^ shows the extinction curve for silicate and 
graphite with the MRN size distribution over amm o, < 

O-max for CLmin — 0.005, Clmax — 0-1 and d^nin — 0.1, dmax — 

0.25. These curves use publicly available extinction data 
calculated using the method of Laor & Draine (1993). 
With the very small grains gone, the graphite absorption 
curve becomes quite flat out to A ~ 1 /im. Figures and ^ 
show the extinction, reddening and mass fraction (relative 
to the full MRN distribution) of dust distributions with 
various value of amm- Curves are given both for (rest- 
frame) E{B — V)/V reddening concentrated at one red- 
shift, and for a cosmological dust distribution (as described 
in §p[). These show that even in the (more reddening) in- 
tegrated extinction, for Qmin ^ 0.06 /im, graphite grains 
give very little {B — V)/V reddening. Sihcate grains do 
not become grey for Omm <, 0.2, but the combined sili- 
cate-|-graphite reddening falls by 50% for amin ^ 0.09 ^m. 
Moreover, this large change in the reddening behavior of 
the dust does not require a large change in the mass: these 
'grey' dust distributions contain 40 — 55% of the mass of 
the MRN distribution. 

6.1. Other Dust Models 

The above conclusions, based on the assumption that 
dust is characterized by the DL model, may not hold for 
other dust models. Mathis & Whiffen (1989) have pro- 
posed that galactic grains are composites of very small 
(a <^ 0.005 ^m) silicate, graphite and amorphous carbon 
particles. These composite grains have a filling factor of 

0.2—1 and corresponding maximal size 0.9 — 0.23 /im. 
Sputtering would be effective at destroying all sizes of low 
filling-factor composite grains, since both gas drag (slow- 
ing the grains) and sputtering would be much more effec- 
tive than in comparably sizes solid spheres. Also, sput- 
tering might tend to 'cleave' large, filamentary grains into 
smaller ones. Large filling-factor particles in this model 



would be much like the Draine & Lee model, although the 
optical properties of the composite materials would differ 
from those of pure graphite or silicate. 

Several core-mantle grain models have also been pro- 
posed; see e.g., Duley, Jones & Williams (1989) and Li 
& Greenberg (1997). The latter model assumes a three- 
component model: large silicate core-organic refractory 
mantle dust, very small carbonaceous grains, and pol- 
yaromatic hydrocarbons (PAHs). The latter two com- 
ponents would presumably be destroyed as dust leaves 
the galaxy, leaving the large core/mantle grains. Li & 
Greenberg take the size distribution of these grains as 
Gaussian, strongly dominated by ~ 0.1 //m grains, with 
parameters chosen to fit the observed extinction curve. 
Such a distribution would be insignificantly affected by 
removal of the small or large-size portions, so intergalac- 
tic dust would have properties exemplified by the large 
core/mantle grains (these grains will redden less than the 
full three-component model, but only slightly). On the 
other hand, it seems that there are good reasons to expect 
a power-law grain size distribution (Biermann & Harwit 
1979; Mathis & Whiffen 1989). It would be interesting to 
investigate whether the model of Li & Greenberg could ac- 
commodate a power law distribution (as they assume for 
the PAHs and the very small grains) . The Duley, Jones & 
Williams (1989) model assumes a bimodal grain-size distri- 
bution: small silicate core/graphite coated grains provide 
UV extinction and the 0.2175 /im bump, whereas an MRN 
distribution of cylindrical silicate grains provides extinc- 
tion in the IR, with grey extinction in UV and optical. In 
this model, intergalactic dust (composed of the large sili- 
cate grains) would be significantly more grey than galactic 
dust, assuming that it can escape. 

Fractal grains (e.g. Wright 1987) and needles (e.g. A99 
and references therein) provide another possible dust com- 
ponent. Needles and platelets have been observed in cap- 
tured dust (Bradley, Brownlee & Veblen 1983), and might 
explain 'very cold' dust in the ISM (Reach et al. 1995). 
As argued in A99 these grains redden very little (especially 
if graphitic) and absorb with high efficiency, hence would 
be preferentially ejected by radiation pressure. Along the 
same lines, DL grains must be at least somewhat elongated 
in order to correctly predict polarization. Elongated grains 
are somewhat more grey than spherical grains of the same 
mass, giving some additional support to the general as- 
sumption that there is a significant grey sub-component 
to interstellar dust. 

While a different grain model might predict a different 
effect of destroying small grains, it is also true that any vi- 
able grain model must be capable of accommodating val- 
ues of Ry 6, since such values are in fact observed. 
Large grains seem to be a necessary component of grain 
models which match the observed extinction laws (Kim et 
al. 1994; Zubko, Krelowski & Wegner 1996, 1998), and 
a greater fraction of these large grains in some regions is 
probably responsible the high observed values of Ry in 
those regions. It is, then, unlikely that the destruction of 
very small grains will make any model more reddening, so 
the assumption of the DL model seems at least qualita- 
tively safe. The MRN grain-size distribution is probably 
safe for the same reasons, and very likely even conserva- 
tive, in the sense that inversions of dust opacity curves into 
grain-size distributions tend to lead to more large grains 
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than MRN would predict (Kim et al. 1994), and the grain- 
size distribution gleaned from observations by the Ulysses 
and Galileo satellites (Frisch et al. 1999) shows many large 
grains up to 1 /im or more in radius. 

7. THE SUPERNOVA RESULTS 

Dust of the DL model with amin ~ 0.1 would corre- 
spond to a dust survival fraction fescfigm ~ 0.4. Using 
equation |, this gives l^XTt ~ (1-5 - 4.5) x 10"^ This 
amount of dust would be quite important cosmologically. 
Measurements of the redshift-magnitude relation of type 
la supernovae (P99; Riess et al. 1998) show statistically 
significant progressive dimming of supernovae which has 
been interpreted as evidence for acceleration in the cosmic 
expansion. This section discusses grey intergalactic dust 
(as specified in f|^) in the context of the these observations. 
Only the DL grain model is considered here. 

Both supernova groups find that after calibration using 
a low-2: sample, the supernova at 2: ~ 0.5 have magni- 
tudes indicative of acceleration in the cosmic expansion. 
The best fit (for a flat cosmology with cosmological con- 
stant) of P99 IS n = 0.28, f^A = 0.72; the results of 
R99 are similar. The necessary extinction to account for 
these results in an = 0.2 open universe (see Fig. ||) is 
Av{z = 0.5) « 0.15 - 0.2 mag.^ An = 1 universe re- 
quires Av{z — 0.5) « 0.4 mag. 

Riess et al. argue that grey extinction would cause too 
much dispersion in the supernova magnitudes to be com- 
patible with their observations if the dust is confined to 
spiral galaxies. Perlniutter et al. derive from their data 
an intrinsic dispersion A at 2; '--^ 0.5 almost identical to 
that at z - 0.05: 

A(z - 0.05) = 0.154±0.04, A(z - 0.5) = 0.157± 0.025. 

This suggests that the processes dominating the intrinsic 
dispersion do not change significantly in magnitude from 
low to high redshift. However, note that - assuming that 
the errors as well as the dispersions add in quadrature - 
the amount of additional dispersion Aadd at high z for- 
mally allowed within the stated errors of is Aadd < 0.13 
mag. This does not include any systematic errors in the 
estimation of the intrinsic dispersion. 

P99 also investigate the mean color difference between 
the high- and low-z samples, finding {E{B — V))zr^o.05 = 
0.033 ±0.014 and {E{B - F ))^^o.5 = 0.035 ± 0.022. Again, 
this suggests that a systematic effect (in color) is not large, 
but nevertheless the errors allow a color difference of 

{EiB - V)),^o.5 - {E{B - V)),^o.05 ^ 0.03 mag. 

In addition, this comparison is subject to a systematic un- 
certainty of ~ 0.03 mag resulting from the conversion of 
(observed) R and / magnitudes into rest-frame B and V 
magnitudes. 

To place tighter constraints on systematic reddening, 
Perlmutter et al. construct an artificially blue subsam- 
ple of the high-z points which is unlikely to be redder (in 
the mean) than the low-z sample. The change in fitting 
that this elimination produces then gives an indication of 
systematic extinction by reddening dust. 

Perlmutter et al. use this method to place a strong con- 
straint of SAv ^ 0.025 mag on the effect of any extinction 

^Ay 0.2 mag accounts for all of the effect, but Ay 0.15 mag 



which (a) exists at high z but not at low z, (b) domi- 
nates the dispersion of both the color and extinction, (c) 
has a reddening-extinction relation Ry up to twice that of 
the Galaxy and (d) occurs in a flat universe. Assumption 
(b), unstated in P99, is crucial but seems unfounded. The 
limit on a systematic increase in dispersion indicates that 
systematic extinction must be a sub-dominant component 
of the total computed 'intrinsic' dispersion in magnitude; 
this holds also for dispersion in color. In this case, remov- 
ing the reddest supernovae will not preferentially remove 
more obscured supernovae, even if dust accounts for the 
whole effect at high z. Thus the stated (more stringent) 
limits on systematic reddening do not apply, as long as the 
dispersion in brightness and/or color is dominated by fac- 
tors other than extinction. Furthermore, if the assumption 
of flatness is dropped, the elimination of the seven reddest 
supernovae actually changes the fit considerably, in the di- 
rection of an open universe (P98, Table 3, Figure 5c). The 
shift corresponds to SAy ~ 0.07 mag at z=0.5. 

The required extinction and the limits on reddening 
and dispersion can now be compared to that expected 
from intergalactic dust. Figure shows the reddening 
Rv for graphite and silicate dust of the DL model, as- 
suming an MRN distribution over Qmin ^ 1 < 0.25 /im. 
The extinction in Figure ^ is integrated to z = 0.5 for 
an Q — 0.2 universe, assuming a constant comoving dust 
density of = 10~^ in each component. The results 

show that, for example, a distribution with amin = 0.1 /im, 
^liusti^ = 0.5) « 4 X 10~^ (total) and equal mass density 
of silicate and graphite grains (approximately the ratio 
derived by Draine & Lee) provides sufficient extinction to 
account for the type la supernova results. The induced 
reddening is 0.025 mag, comparable to the allowed red- 
dening due to either random or systematic errors. Most of 
this reddening is provided by the silicate grains, so the 1:1 
ratio is conservative; the real ratio should be biased toward 
the less efficiently destroyed (Draine & Salpeter 1979a) and 
possibly more efficiently ejected (assuming radiation pres- 
sure expulsion) graphite grains. The large grains contain 
^ 40% of the full MRN distribution. Larger values of amin 
provide less reddening but values of amin ^ 0.15 proba- 
bly contain too little mass to be viable in explaining the 
supernova data. Graphite grains alone (if silicate grains 
were preferentially destroyed) with amin 0.06 (giving 

fescfigm <, 0.6) and ~ 3 X 10~^ would produce sim- 

ilar effects. 

The amount of dispersion induced by the dust is very 
important but can be estimated only roughly. Assum- 
ing that the dust is uniformly distributed in randomly 
placed spheres of radius R with number density n, the 
dispersion A is given approximately by A/Av{z = 0.5) « 
N~^/^, where N is the number of spheres intersected 
by a typical path, and can be written TV ~ mrR^D, 
where D « 2400ft.^g^ Mpc is the distance to z = 0.5 
in an Q = 0.2 universe. Now consider galaxies with 
n = OmShl^il + z)3Mpc-3 (Lin et al. 1996). For 
^ 1 this implies R ^ 70h^l[{l + z)/1.5]-^^^ kpc. 

Escape velocities from spirals are ^ 250kms~^, so the 
dispersion in integrated optical depth due to dust ejected 



;s an n = 0.2 open universe within the stated I — a contour. 
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by radiation pressure or winds and traveling away from 
the disk for time r^g^ is 

Figure ^ shows the (small) difference in optical depth be- 
tween the total dust distribution and the dust which has 
existed for > 200 Myr and > 1 Gyr, demonstrating that 
the dispersion induced by grey dust created but not yet 
sufRciently dispersed would be small. Large-scale correla- 
tions between galaxies are unlikely to be important in this 
analysis. The dispersion in density on 8 Mpc scales is ^ 1, 
but ~ 300 such do mains lie between here and z ~ 0.5, 
leading to A - 0.2/V300 = 0.01. 

The best way to estimate dispersion (currently under- 
way) would probably be to measure the optical depth 
through random lines of sight piercing a high-resolution 
cosmology simulation which either tracks metals (Cen & 
□striker 1999b; Gnedin 1998) or allows some prescription 
relating gas density to dust density in the IGM. Short of 
this, I note that Cen & Ostriker's simulations show that 
5 ^ 100 is characteristic of the bulk of the metal-rich gas. 
This corresponds to i? ^ 670ft.g5^(l + z)~'^ kpc for the 
spheres considered above. Fairly uniform regions of this 
size and number density would give little dispersion. 

If i? 70 kpc, most extinction takes place in a small 
number of clumps, in particular in the halo of the super- 
nova host galaxy. For dust of density ^dust = X ^ 10~^ 
uniformly distributed in radius i? = ^ x 100 kpc halos of 
galaxies with n = 0.008/ig5(l + zY Mpc~^, the extinction 
to the galaxy center is 

= T^n 

= »-"^( 5xlOWg-0 '^'-'"^"'' 

which gives the required extinction for ^ « 1 and % « 4. 
As long as ^ ^ 0.1 there will not be large dispersion due 
to different galactic radii in the supernovae, but inhomo- 
geneities in the dust distribution in the halo could be im- 
portant. The supernova results could could be explained 
by such halos with a smaller ^dust , but this would require 
that the halos are somewhat larger at low z, and that 
all halos at z ~ 0.5 have similar column density through 
them. This scenario does not seem to be as natural an 
explanation as a more uniform dust distribution, but it 
remains a possibility. 

Leaving aside these uncertainties, the essential result of 
my calculation is that a truncated-MRN distribution of 
Draine & Lee dust with amin 0.1, uniformly distributed 

and with ^^dust{z — 0.5) ~ 4 x 10^^ can account for the su- 
pernova dimming in an f2 = 0.2 universe without excessive 
reddening. Whether the induced dispersion in extinction is 
too large depends crucially on the (uncertain) distance to 
which ejected dust can escape from the galaxies in which it 
forms, and on the dumpiness of the resulting distribution. 

^Intergalactic dust in the calculations of Loeb & Haiman (1997) 
UV/optical background lower than recent detections indicate. 



8. COSMIC BACKGROUNDS 

An intergalactic dust distribution of the magnitude re- 
quired to explain the supernova results would have other 
cosmological implications. For instance, any intergalactic 
dust component will absorb energy from the optical/UV 
background and re-emit the energy in the FIR/microwave. 
Calculation of the evolution of the cosmic mean density 
field shows (Aguirre & Haiman, in preparation) that the 
dust considered in this paper will not lead to measurable 
CMB spectral distortions,^ but instead adds (significantly) 
to the GIB. The FIRAS distortion limits can, however, 
limit dust types with high FIR emissivity, since these have 
lower equilibrium temperatures. 

Note also that while the SFR estimates correct for dust 
extinction (e.g. Madau, Pozetti & Dickinson 1998; Pettini 
1999 and references therein), these corrections would not 
account for intergalactic dust. Intergalactic absorption of 
~ 0.1 — 1 mag at z ~ 0.5 — 5 (see fig. ||) would imply an 
SFR - and hence metal density - a factor of ^ 2 higher 
than that given in 

9. TESTING FOR INTERGALACTIC DUST 

Future observations of supernovae can investigate the 
importance of intergalactic dust in two ways. First, ac- 
curate high-z observations in rest-frame i?-band or longer 
wavelengths should reveal the dust with properties of the 
model developed here. As shown in fig. |[ the E{B — R) / B 
reddening to z = 0.5 is ~ 0.25 for both pure graphite 
dust (with a„iin — 0.06 /im) or for silicate-|-graphite (with 
o-min = 0.1 /im). Therefore high-z supernova should have 
E{B — R) values 0.05 mag higher than the low-z sam- 
ple. The effect \x\ B — I would be even stronger. Non- 
spherical grains (of which the needles of A99 are an exam- 
ple) could provide grey opacity into i?-band and beyond, 
but are constrained limits on far- infrared/microwave emis- 
sion (Aguirre & Haiman, in preparation). 

Figure || includes the deviation of the fit of R99 from the 
SI = 0.2 model (dashed line). The difference between this 
line and the dust optical depth then indicates the deviation 
at high z of the dust model from the cosmological constant 
(f^A ~ 0.72) model; at z = 1.5 the difference is ~ 0.2 — 0.3 
mag. This possible difference should be testable once a 
significant number of z > 1 supernovae are measured. A 
large number of supernovae is important both because the 
dispersion is comparable to the effect being measured, and 
because it is crucial to be sure that the z > 1 sample is 
statistically complete ~ a worry that is less significant at 
z 0.5 since the claimed effect is of dimming rather than 
brightening. 

An intergalactic dust distribution would also dim other 
distant objects. Dust has been proposed several times as 
an explanation for the dropoff in quasar number counts for 
z ^ 3 (e.g. Ostriker & Heisler 1984; Wright 1986). While 
this point is still controversial, observation of the dropoff 
in radio-selected quasars (Shaver et al. 1996) shows this 
explanation to be unlikely. The models of this paper would 
not predict the dropoff (unless substantial population HI 
dust exists also) since Av{z) <^ 1 for all z. Galaxy counts, 
in contrast to quasar studies, show an excess in i3-band 
counts at high z as compared to some evolution models 

Ferrara et al. (1999) has low temperature, but both papers assume a 
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(see Shimasaku & Fukugita 1998 for a summary). For 
B-band observations probing z ~ 0.5 ± 0.3, this 'excess' 
amounts to a deviation of up to '--^ 1 mag from the pre- 
dicted curve. The dust extinction of <J 0.4 mag for z ^ 1 
would enlarge this discrepancy, but not greatly. This holds 
for / and if -band surveys as well. Because of uncertain- 
ties in the models at high z, Shimasaku & Fukugita are 
reluctant to draw conclusions about cosmological param- 
eters from the galaxy counts. But with rapid progress in 
the field this may soon become a useful test of the dust 
model. 

Several claims of dust detection in clusters have been 
made, but remain controversial. The model of A99 pre- 
dicted that cluster dust should be more grey than galactic 
dust; the arguments of this paper reinforce this predic- 
tion. The intracluster gas destroys dust efficiently; since 
sputtering destroys small grains more efficiently than large 
grains, it is likely that whatever grains survive in clusters 
are large, and thus supply grey opacity (this may help 
explain the controversy surrounding the existence on in- 
tracluster dust). Determination of the extinction curve of 
dust in clusters (if dust indeed exists) would be an impor- 
tant test of the key idea of the intergalactic dust model, 
although grey dust in clusters would not necessarily imply 
that dust in the diffuse IGM is also grey. 

Dust confined to galactic halos might also be detectable 
by its IR emission or by studies of objects seen through the 
halo. Zaritsky (1994) has presented preliminary evidence 
for halo dust in NGC 2835 and NGC 3521 at ~ 60kpc 
using B — I reddening; if confirmed this would be an im- 
portant dust component, as the inferred mass is large. 

10. CONCLUSIONS 

As discussed in §^ fairly strong arguments suggest that 
the universe currently has ilz ~ (1.5 — 3) x 10"'* in metals. 
This metal cannot all be contained in the stars and gas in 
known galaxies, so unless metals are well hidden in unob- 
served galaxies, a metal density $7^™ « (0.75—2.25) x lO"'* 
should exist in the more uniform IGM or in extended ha- 
los. Moreover, most of this metal was probably in place 
by z ~ 0.5. 

Some fraction d,n of this metal must be in dust. There 
are several good reasons to expect that the grain size dis- 
tribution of this dust should be different than for dust in 
galaxies: (1) Ejection by radiation pressure favors higher 
opacity, less reddening grains (A99). Within the DL 
model, this means large grains are preferentially ejected. 
(2) As dust leaves galaxies (where it is almost certainly cre- 
ated) , small grains are preferentially destroyed by sputter- 
ing. This preferential destruction occurs also in the IGM, 
but the efficiency is rather uncertain. (3) Small grains are 
generally assumed to form from the shattering of larger 
grains. This shattering will not occur in the IGM due 
to the low densities, nor can small grains grow from the 
vapor. 

The principal assumption of this paper (supported by 
what detailed calculations are available) is that very small 
dust grains leaving the galaxy - comprising a fraction 
(1 — /) of the total dust density - are removed from the 
grain-size distribution, while large grains are not. Dust 
as modeled by DL, with grains of radius (^ 0.1 /im) re- 
moved (/ « 0.4), reddens very little yet has a visual opac- 



ity k w 5 x lO'^cm^g^^. Uniformly distributed dust of 
constant comoving density in an Q = 0.2 open universe 
provides an extinction to z = 0.5 of 

Av 0.15/165 I "tn" 5 I 104"' 2 1 ) ' (9) 

\ 4 X 10"^ / \5 X 10 cm"' g^^ / 

hence it can account for the dimming of type la supernovae 
at z ~ 0.5 in a way fully consistent with observations. 
The fact that the expected dust density of 0.5/f2^™ ~ 
(0.25 - 2.25) X 10"'' X (0.5) x (0.4) = (1.5 - 4.5) x 10"^ 
is so close to the density required is very interesting. Ac- 
cording to these estimates (see also §0 and §7) it is pos- 
sible, but rather unlikely, that there is sufficient dust 
i^liusti^ — 0-5) ~ 9 X 10"^) to allow compatibility with a 
closed, matter-dominated universe. 

The arguments of this paper show that there may be an 
intergalactic dust component which is important yet has 
evaded earlier attempts at its discovery; from this point of 
view the supernova observations are telling us not about 
cosmic acceleration, but about cosmic opacity. If the argu- 
ments of this paper are substantially correct, an intergalac- 
tic dust distribution of some magnitude is inevitable, so it 
is worth pointing out the chief caveats (which maintain 
rough consistency with other observations) under which 
the obscuration would not be sufficient to to account for 
the supernova dimming: 

• Both the SFR argument and the cluster enrich- 
ment argument may predict cosmic metallicity sev- 
eral times the true value. In this case most metal 
could be locked in galaxies. Taking this view, how- 
ever, probably also requires one to assume that clus- 
ter galaxies both have a different IMF and eject met- 
als more efficiently than field galaxies do. 

• A large population of unobserved but compact ob- 
jects such a low surface brightness galaxies could 
contain a large fraction of the cosmic metallicity. 
This would, however, also fail to explain the clus- 
ter observations. 

• The Draine & Lee model may not accurately char- 
acterize dust, so the effect of small grain destruction 
on the opacity curve may be significantly different 
than assumed here. But the alternative grain model 
would then probably have to exclude the grey sub- 
component which most current grain models include, 
yet still explain the values Ry ^ 6 observed in some 
regions. 

• Dust of all sizes may be efficiently destroyed in the 
diffuse IGM or as it leaves galaxies. While total grain 
destruction does not seem to be quantitatively sup- 
ported in the DL model, it is certainly plausible, and 
would probably be the case if grains are very fluffy, 
high filling-factor composites. 

It is unlikely that the issue can be cleanly decided on the 
basis of the above points or by the current supernova data 
(i.e. limits on dispersion or reddening); on the other hand, 
future supernova results can rule decisively. Statistically 
robust deviation of the magnitude-redshift curve of z > 1 
supernovae from the dust prediction would argue strongly 
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for the interpretation of the dimming as cosmic accelera- 
tion and for the relative unimportance of grey intergalac- 
tic dust. Clear evidence of, for example, systematic (rest 
frame) B — R reddening would argue strongly for dust. 

If future observations show that intergalactic dust is in- 
deed important, observations of type la supernovae will be 
seen to have not only determined the deceleration param- 
eter (once dust is accounted for), but to have discovered 
a component of the universe that will have important im- 
plications for many future observations at high redshift. 
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INTERGALACTIC DUST AND TYPE lA SUPERNOVAE 




Fig. 1. — Comoving star formation rate (thick; left axis) and integrated (from z = 10) cosmic metal density (thin; right axis), from 
Madau (1999; solid) and Steidel et al. (1999; dashed). Adjusted for Hq = Gbkms'^ Mpc'^, fl = 0.2. 
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Fig. 2. — Opacity (absorption+scattering) of Draine & Lee model with truncated MRN distribution amin — 0.005 /im, amax = 0.1 for 
graphite (sohd, thin) and silicate (dashed, thin), and for amin = 0.1 /im, amax = 0.25 for graphite (solid, thick) and silicate (dashed, thick). 
Vertical dotted lines indicate B and V-band centers. 
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Fig. ,3. — Top: Extinction to z = 0.5 in f2 = 0.2 universe for {z = 0.5) = 10 in graphite (solid), silicate (dashed) and graphite+silicate 

(dot-dashed). Bottom: Fraction of full MRN distribution contained the distribution with amin < a < 0.25 jum. 
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Fig. 4. — Top: E(B — V)/V = Ry^ reddening for graphite (solid) and silicate-fgraphite (dashed), for constant comoving dust density 
integrated to z = 0.5 with Q = 0.2; E{B — R)/B reddening for graphite (dot-dashed) and graphite+silicate (triple-dot-dashed). The filters 
refer to rest frame colors at z = 0.5. Bottom: Same, but for dust at a single redshift with colors referring to that frame. 
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Fig. 5. — Deviations from an open universe: Optical depth to redshift ^ in (2 = 0.5 rest-frame) V-band (solid, dark). Deviation of the 
P99 fit = 0.72, f2 = 0.28 from a O = 0.2 open cosmology (dashed). Solid, light lines indicated optical depth due to dust formed dt = 200 
Myr and dt = 1000 Myr earlier than indicated by dark solid line. The dotted line signifies the 'intrinsic' dispersion in observed supernova 
magnitudes about the fit. The assumed dust density is proportional to the integrated (Madau 1999) SFR (see figure hi), normalized to 
Q'/'"(2 = 0)=4.5xlO-^ 

dust ^ ' 



